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INTRODUCTION 


The  University  of  Massachusetts  Lowell  Center  for  Atmospheric  Research  (UMLCAR) 
continues  to  work  towards  a  better  vmderstanding  of  the  physics  of  the  upper  atmosphere  that  can 
lead  to  reliable  transionospheric  and  subionospheric  communications,  meeting  the  needs  of  the 
United  States  Air  Force.  With  the  arrival  of  networking  and  super  communication  highways, 
space  science  has  been  propelled  into  real-time  forecasting.  New  capabilities  will  greatly 
enhance  the  ability  of  the  United  States  Air  Force  to  determine  when  communication  systems 
may  be  affected  by  ionospheric  and/or  solar-terrestrial  phenomena.  These  capabilities,  being 
pursued  at  UMLCAR  in  conjunction  with  other  institutions  as  part  of  the  “Space  Weather 
Forecasting”  initiative,  are  essential  for  any  program  designed  to  use  the  space-earth 
environment  for  reliable  future  communication  requirements. 

The  task  of  specifying  the  ionosphere  has  been  subdivided  into  several  areas.  These 
include: 


■  Structure  of  Equatorial  Depletion  Bands: 

The  plan  is  to  study  the  morphology  of  equatorial  depletions  and  to  determine  the 
physical  cause  for  their  development  and  behavior.  Particular  emphasis  will  be  placed 
on  the  use  of  ionosonde  and  other  high  frequency  (HF)  measuring  techniques 
involving  both  the  experimental  and  theoretical  study  of  these  phenomena.  The 
experimental  part  consists  of  measurements  made  with  a  variety  of  instruments 
during  the  campaigns  planned  by  Air  Force  Research  Laboratoiy  (AFRL)  scientists  as 
well  as  using  the  routine  observations  at  established  equatorial  sites  such  as  at 
Jicamarca,  Peru. 

■  Polar  Cap  Monitoring: 

Ionospheric  drift  measurements  from  high  latitude  Digital  Ionospheric  Sounding 
System  (DISS)  stations  have  shown  some  correlation  between  the  sign  of  Bz  and  the 
plasma  drift  velocity.  Ionospheric  drift  data  from  digital  ionosondes  at  Qaanaaq, 
Greenland;  Sondrestrom^ord,  Greenland;  Svalbard,  Norway;  and  Goose  Bay, 
Labrador,  are  used  to  refine  the  process  of  determining  the  polarity  of  Bz  from 
ground-based  ionospheric  drift  measurements.  These  results  are  used  to  develop  the 
polar  cap  convection  patterns. 

■  Ionospheric  Modification  Diagnostics: 

Under  this  task,  UMLCAR  will  be  responsible  for  deploying  the  ionospheric 
diagnostic  equipment,  conducting  measurements  and  providing  quick-look  data  for 
the  conduct  of  ionospheric  heating  experiments  and  analysis.  UMLCAR  also  provides 
geophysical  interpretation  of  the  data. 

■  Field  Campaigns: 

Under  this  task  UMLCAR  participates  in  measurements  at  remote  field  locations. 
UMLCAR  provides  support  to  AFRL  in  making  joint  measurements  with  other 
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sensors  using  tiie  DPS  systems.  As  part  of  this  task,  UMLCAR  will  construct  one 
four-receiver  Digisonde  Portable  Soimder  (DPS-4)  system  with  four  receive  antennas. 

■  Development  of  Digisonde  Database: 

UMLCAR  collects  and  provides  quality  control  on  electron  density  profile  data  for 
selected  periods  and  locations.  These  refined  data  are  used  to  compare  with  the  Air 
Force  Research  Laboratoiy’s  Parametrized  Ionospheric  Model  (PIM).  Also, 
UMLCAR  is  developing  a  standard  format  for  these  analyzed  Digisonde  data 
and  archives  these  data  in  the  Digisonde  database. 

■  Digisonde  Network  Support: 

UMLCAR  provides  upgrading  of  DISS,  AN/FMQ-12  and  network 
support,  particularly  for  building  the  ionospheric  model  database. 

■  Low  Power  D^tal  Portable  Sounder  (LPD): 

UMLCAR  is  developing  a  prototype  of  a  Low  Power  Digital  lonosonde 
(less  than  100  watts)  needed  for  routine  ionospheric  measurements. 

■  Topside  Automated  Sounder,  TOPAS: 

Under  this  task,  UMLCAR  is  designing  and  building  a  compact  and 
efficient  prototype  ionosonde  for  satellite  based  observations. 

■  Propagation  Analysis  (Project  SEND): 

Using  analytic  ray  tracing  techniques,  UMLCAR  was  tasked  to  develop 
a  program  for  rapid  propagation  analysis  with  expected  signal-to-noise 
as  an  indicator  to  the  USAF  user  of  the  expected  performance  of  his  HF 
system. 

2.  RESEARCH 

2.1  Structure  of  the  Equatorial  Ionosphere 
2.1.1  Equatorial  Depletion  Analysis 

Before  attacking  the  dynamics  of  the  depletion  regions  observed  at  Agua  Verde,  Chile,  it 
was  decided  to  investigate  the  chemistry  and  dynamics  of  the  background  F-layer  as  it  affects  the 
overall  630.0  nm  emission  intensity.  Following  our  earlier  work  involving  the  chemistry  and 
vertical  motions  within  the  F-layer  that  separated  the  contribution  of  each  to  the  observed 
Doppler  shift,  it  was  decided  to  apply  this  technique  to  determine  the  integrated  airglow  emission 
intensity. 
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A  rather  simple  model  was  developed,  limited  at  this  stage  to  post  F-region  sunset  times 
when  it  could  be  assumed  that  there  was  no  solar  production  of  ionization.  Under  these 
conditions  in  the  F-region,  the  chemistry  becomes  rather  simple  and  only  an  electron  loss 
mechanism  was  considered  that  involves,  jBrst,  charge  exchange  between  and  O2,  that 
followed  by  a  dissociative  recombination  reaction  which  removes  the  electrons  and  produces 
excited  atomic  oxygen.  The  O*  in  a  forbidden  state,  after  some  time,  emits  a  “red”  photon  and 
the  atomic  oxygen  returns  to  the  ground  state.  It  is  these  photons  that  constitute  the  observed 
630.0  nm  airglow. 

In  addition  to  the  loss  of  electrons  by  this  recombination  process,  the  redistribution  of 
electrons  by  vertical  transport  was  considered.  Both  upward  and  downward  transport  were 
considered.  The  electron  density  profiles  as  a  fimction  of  time  were  calculated  using  the 
continuity  equation  including  the  effects  of  vertical  transport.  For  this  problem,  changes  in  the 
electron  density  in  the  250  to  300  km  altitude  regimes  were  critical  as  the  somce  region  of  the 
airglow.  Typically,  this  altitude  range,  in  the  early  evening,  is  at  or  below  the  F-layer  peak. 
Assuming  that  the  reaction  rates  for  the  chemistry  part  of  this  problem  are  relatively  well  known 
(these  simulations  were  run  using  the  reaction  rates  taken  from  the  rather  complete  airglow 
emission  program  developed  by  Stan  Solomon  at  the  University  of  Colorado),  the  measured 
electron  density  profiles  were  fit  by  varying  the  vertical  transport  velocity  as  a  fimction  of  time 
as  the  input  to  the  continuity  (differential)  equation.  The  assumption  was  made  that  the  transport 
velocity  is  constant  over  the  altitude  range  for  this  simulation,  i.e.,  200  to  600  km.  The  altitude 
range  for  integration  of  the  differential  equation  was  made  larger  than  the  region  of  interest, 
thereby  moving  the  boimdary  conditions  out  of  the  critical  region. 

For  October  1,  2,  and  3,  1994,  the  measured  plasma  frequency  contoms  at  Agua  Verde  as 
a  function  of  time  are  shown  below. 


Agua  Verde,  Chile  -  Oct  2, 1994  Agua  Verde,  Chile  -  Oct  3,  1994 
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Tlme(minutBSx5) 

Agua  Verde,  Chile  -  Oct  3, 1994 

Figure  1.  Plasma  frequency  contours  for  October  1, 2,  and  3,  1994. 

The  contours  are  in  0.5  MHz  steps.  Time  axis  begins  at  00  UT. 

As  stated  above,  die  differential  equation  requires  a  vertical  velocity  time  profile 
(assumed  uniform  widi  altitude)  and  after  several  trials  the  profile  in  Figure  2  was  used  to 
generate  a  series  of  electron  density  height  profiles  (die  step  intervals  for  the  solution  of  the 
differential  equation  were  1  s  and  5  ^). 


Time  (seconds  after  0030  UT) 


Figure  2.  Assumed  vertical  velocity  time  profile  for  October  2, 1994,  at  Agua  Verde. 
The  time  scale  begins  at  0030  UT.  (Downward  is  negative) 


4 


z 


Figure  3.  Calculated  electron  density  profiles  at  0, 40, 60  and  90  minutes  after  0030  UT. 

A  downward  velocity  of  10  m/s  was  assumed  to  persist  for  30  minutes,  followed  by  a  30- 
minute  transition  to  an  upward  motion  of  7  m/s.  The  zero  vertical  velocity  was  set  at  t  =  3000s, 
i.e.,  at  0120  UT  for  these  particular  data.  The  results  of  this  computation  (Figure  3)  show  a 
selected  sequence  of  electron  density  profiles  at  t  =  0,  40,  60  and  90  minutes  after  die  start  at 
0030  UT. 

The  volumetric  emission  rate  (photons/cm^/s)  was  calculated  using  the  rate  of  loss  of 
electrons  at  each  altitude.  The  results  for  October  2, 1994,  at  0030  UT  are  shown  in  Figure  4.  In 
this  case,  the  maximum  emission  rate  of  just  below  100  (photons/cmVs)  occurs  at  about  240  km 
with  a  relatively  rapid  fall  off  at  lower  altitudes  and  a  more  gradu^  decrease  above.  The 
calculated  emission  rate  is  down  by  almost  a  factor  of  ten  at  ±  50  km  fi-om  the  peak  altitude. 
Finally,  for  each  time  (one-second  st^s)  the  volumetric  emission  rate  is  integrated  over  the 
altitude  range  of 200  to  600  km  to  produce  calculated  emission  intensity  for  the  simulation  of  the 
measured  optical  airglow  on  the  ground.  These  results  are  shown  in  Figure  5. 


Height  (km) 

Figure  4.  Calculated  volumetric  emissions  rate  (630.0  nm)  at  0030  UT. 
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Time  (minutes  after  0030  UT) 


Figure  5.  Calculated  emission  intensity  of  the  630.0  nm  airglow  for  October  2, 1994,  at  Agua  Verde 

based  on  simulated  electron  density  profiles. 

Here,  the  emission  intensity  begins  at  570  R  and  decreases  rather  steadily  over  90 
minutes  to  about  180  R.  Comparing  this  curve  with  the  measured  emission  intensity,  fi-om  30 
minutes  to  90  minutes  after  00  UT,  the  general  behavior  is  very  similar.  The  obvious  and 
disturbing  part  of  this  comparison  was  the  discrepancy  in  the  absolute  value  of  the  emission 
intensity  measured  in  Rayleighs,  differing  by  more  fiian  a  factor  of  5. 

To  investigate  the  problem  of  the  discrepancy  between  the  calculated  and  the  measured 
emission  intensity,  it  was  decided  to  consider  several  different  approaches. 

First,  by  rerunning  the  program  under  different  conditions,  the  theoretical  analysis 
developed  at  UMLCAR  was  reexamined.  At  that  point,  only  the  consistency  of  the  several 
solutions  was  used  to  indicate  that  the  program  was  operating  properly.  Given  these  reasonable 
results  and  after  discussions  with  E.  Weber,  it  was  decided  to  take  Ae  following  two  directions. 

1.  Stan  Solomon  was  contacted  at  the  University  of  Colorado.  He  has  developed  an 
extensive  code  to  calculate  the  intensity  of  630.0  nm  emissions  in  the  F-region. 
His  code  is  considerably  more  detailed  than  the  one  developed  here.  A  meeting 
was  held  with  Solomon  after  an  AFRL  seminar  and  it  was  agreed  that  an  electron 
density  profile  firom  Agua  Verde  would  be  sent  to  him  for  analysis  using  his 
emission  code.  It  then  would  be  possible  to  compare  his  results  with  our  results. 
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2.  Meetings  were  held  with  Peter  Ning  and  his  associates  from  KEO  to  discuss  the 
calibration  of  the  HAARP  scanning  photometer.  It  was  agreed  that  they  would 
look  into  possible  sources  for  the  discrepancy  between  the  CAR  simulations  and 
the  “measured”  emission  intensities  that  were  larger,  by  a  factor  of  5  or  6. 

Working  with  Solomon  to  get  the  details  straightened  out,  the  calculated  emission 
intensities,  using  his  elaborate  code,  were  foimd  to  be  within  10%  of  our  calculated  value. 
Almost  simultaneously,  the  report  came  back  from  KEO  that  they  had  made  a  mistake  when  they 
believed  that  they  had  “averaged”  four  adjacent  pixels  to  reduce  the  512  x  512  array  of  CCD 
elements  to  256  x  256.  Actually,  the  HAARP  system  siunmed  the  four  pixels  and  did  not 
average  them. 

Doing  the  averaging  increased  the  “measured”  intensity  by  a  factor  of  four  and  we  feel 
that  the  basic  discrepancy  has  been  resolved  and  have  concluded  that  our  simple  emission  code 
including  vertical  drift  motions  to  follow  the  development  of  the  layer  and  the  associated  airglow 
emission  intensity  after  sunset,  performs  well.  The  remaining  difference  between  the  measured 
and  calculated  intensities,  typically  30  to  50  percent,  needs  to  be  addressed  and  may  involve 
adjusting  the  electron  temperature  or  one  of  the  several  other  factors  in  the  chemistry  of  the 
simulation. 

With  a  simulation  code  that  works  reasonably  well,  it  is  now  possible  to  investigate  the 
dynamics  of  the  depletion  bands,  applying  the  same  code  with  different  vertical  transport 
velocity  profiles  to  determine  the  range  of  speeds  that  lead  to  partially  depleted  electron  density 
bands.  This  research  addresses  the  problem  of  the  residual  emissions  consistently  observed 
inside  the  depletion  bands  at  Agua  Verde. 

After  resolving  the  question  of  the  calibration  of  the  photometer  and  achieving  good 
agreement  between  the  observed  emission  rates  and  calculations  using  the  measured  electron 
density  profiles  and  the  simplified  version  of  the  chemical  and  transport  code,  the  investigation 
of  the  observed  emission  inside  the  depletion  and  the  variations  of  the  background  emission  rate 
spatially  over  the  sky  and  from  day  to  day  was  begun.  The  data  used  included  both  the 
Digisonde  and  the  all-sky  photometer  images  at  630.0  run  and  focused  on  two  problem  areas. 

1.  Understanding  the  differences  between  the  emission  rate  on  October  2  when  no 
depletions  were  present  (the  observed  rate  was  about  2x  the  background  rate 
measured  on  October  1  and  3)  and  on  the  other  nights  when  several  depletions  were 
present  during  the  pre-midnight  period. 

2.  On  both  nights  when  depletions  were  present,  there  appears  to  be  a  depression  of  the 
ionosphere  of  the  order  of  20  to  30  km  before  the  moving  depletion  arrives  at  the  site. 

The  soimder  and  the  all-sky  data  were  reexamined  to  confirm  these  observations.  These 
results  were  presented  at  the  AGU  spring  meeting  in  Boston  and  at  the  International  Union  of 
Radio  Science  (URSI)  meeting  in  Toronto  in  August. 
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2.1.2  Equatorial  Spread-F 


The  paper  “Pacific  Region  Equatorial  Anomaly  Studies  in  Asia  (PREASA)”  written  by 
J.L.  Scali,  K.  Igarashi,  B.W.  Reinisch»  H.  Minakoshi  and  Y.  Masuda,  has  been  submitted  to  the 
Joiunal  of  Geophysical  Research  Space  Physics  (JGR)  for  publication.  In  this  p^er,  ionospheric 
irregularities  associated  with  a  daytime  scintillation  event  observed  in  the  Japan  Asian  sector 
during  the  PREASA-H  campaign  on  23  February  1997  are  described.  The  scintillation-causing 
irregularities  are  identified  as  small-scale  field  aligned  depletions  in  the  lower  F  region.  The 
duration  of  scintillation  is  controlled  by  how  rapidly  these  irregularities  are  filled  in  by  daytime 
ionization  sources.  The  triggering  mechanism  appears  to  be  related  to  the  occurrence  of  spread 
sporadic  E  and  displays  propagation  speeds  from  high  to  low  latitudes  comparable  to  medium- 
scale  gravity  waves. 


2.2  Propagation  Analysis  (Project  SEND) 

At  a  recent  meeting,  the  USAF  established  a  new  initiative  called  Project  SEND.  The 
requirement  was  to  develop  an  HF  propagation  package  that  was  capable  of  modeling  signal, 
noise  and  signal-to-noise  ratio  (SNR)  for  the  coverage  area  of  a  communications  system  at  any 
arbitrary  location  on  the  surface  of  the  earth.  The  approach  was  to  use  the  Air  Force’s  adjustable 
ionospheric  model,  PRISM,  in  connection  with  a  rapid  ray  tracing  routine.  In  the  particular 
theater  of  operation,  where  the  “wanted”  transmitter  is  located,  it  was  assumed  that  a  DISS 
(vertical  sounder)  would  be  available  to  update  PRISM  and  thereby  provide  a  “real-time” 
ionosphere  for  the  ray  tracing. 

The  goal  was  to  compute  the  expected  signal  level  and  SNR  over  a  360  degree  azimuthal 
scan  that  typically  would  be  provided  by  the  wanted  transmitter.  These  computations  were  to  be 
available,  potentially;  to  several  customers  in  near  real-time  (approximately  every  15  minutes). 
It  was  also  necessary  to  consider  the  effect  of  an  “unwanted”  transmitter  operating  at  the  same 
frequency  as  the  wanted  transmitter,  effectively  trying  to  jam  any  receiving  system  operating  at 
the  wanted  frequency.  This  “jammer”  was  also  to  have  an  arbitrary  location  within  the  same 
coverage  of  the  PRISM  model  as  updated  by  DISS.  It  was  decided  that  the  PRISM  area  in  which 
the  wanted  and  unwanted  transmitters  were  located  would  be  a  4000  km  radius  circle  around  the 
DISS  site. 

HASEL,  which  is  the  3-D  numerical  ray  tracing  code  developed  by  Chris  Coleman  at 
DSTO,  Australia,  was  the  first  choice.  Although  this  program  has  great  precision,  it  runs 
relatively  slowly,  particularly  for  this  project  where  it  might  be  required  to  calculate  the  coverage 
for  several  frequencies  for  several  customers.  Each  frequency  required  some  180  azimuthal 
steps  and  up  to  90  elevation  steps  or  approximately  15,000  rays  per  frequency  per  customer. 
HASEL  was  too  slow  and  an  analytic  ray  tracing  code  “SM^T”  developed  at  the  Defense 
Research  Establishment  in  the  UK  was  considered. 

SMART  is  an  analytic  ray  tracing  code  that  assumes  each  altitude  slice  is  quasi-parabolic 
in  shape  and  it  is  then  possible  to  have  a  rapid  calculation  of  the  ray  path  in  that  segment.  Each 
segment,  in  altitude  and  range  was  modeled  in  a  quasi-parabolic  form  and  the  calculated  ray 
paths  were  extended  from  the  transmitter  up  to  the  ionosphere  and  back  to  the  ground.  On  a 
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UNIX-based  SGI  computer,  each  frequency  (15,000  rays)  was  calculated  in  only  a  few  minutes, 
sufficiently  fast  to  meet  the  operational  requirements. 

SMART  had  two  other  advantages.  First,  it  calculated  the  geometrical  spreading  loss  for 
one  and  two  hops  (two  hops  are  required  to  get  the  transmitter  coverage  out  to  near  4000  km); 
and  second,  it  calculated  the  absorption  of  the  radio  wave  through  the  lower  ionosphere:  neither 
of  these  features  was  available  with  HASEL.  SMART  also  included  a  grid  mapping  routine  that 
took  the  radial  array  of  rays  from  the  transmitter  and  generated  a  pixel-by-pixel  map  for  the 
signal  strength  or  the  SNR. 

Tests  were  run  to  compare  the  performance  of  SMART  against  the  more  accurate 
HASEL.  For  the  relatively  slowly  varying  ionospheres  generated  by  PRISM,  the  SMART 
results  were  sufficiently  close  to  HASEL  and  the  listed  advantages  made  the  SMART  code  the 
smart  choice.  On  the  other  hand,  SMART  did  not  have  a  noise  model  and  after  some  research 
the  NOISE.DAT  model  distributed  by  the  International  Telecommunications  Union  was  obtained 
and  incorporated  as  part  of  the  SEND  software  package.  This  computer-based  model  was  the 
updated  version  of  CCIR  312,  which  served  the  radio  propagation  community  for  many  years. 
This  noise  model  gives,  for  any  location,  season  and  time  of  day,  the  atmospheric  noise  level  at 
the  selected  frequency.  It  also  provides  the  statistics  associated  with  the  predicted  noise  level. 
This  noise  code  was  incorporated  into  SMART,  which  then  made  it  possible  to  compute  the  SNR 
at  each  pixel  within  the  coverage  map.  NOISE.DAT  was  also  able  to  estimate  the  galactic  noise 
component,  but  galactic  noise,  typically,  lies  below  atmospheric  noise  for  frequencies  below  30 
MHz.  Finally  NOISE.DAT  estimates  the  man-made  component  of  noise  depending  on  the 
particular  environment.  This  has  to  be  selected  from  four  choices,  beginning  with  the  quietest: 
quiet  rural,  to  rural,  to  residential  and  finally  to  the  highest  noise  level,  industrial.  Unfortunately, 
this  is  a  subjective  choice,  depending  on  the  location  of  the  receive  site.  Without  being  able  to 
specify  a  specific  receive  site,  the  default  manmade  noise  level  was  chosen  to  be  residential  and 
applied  ever5where. 

Finally,  it  was  decided  to  modify  the  absorption  model  in  SMART.  The  developers  of 
SMART  used  the  typical  control  point  approach,  i.e.,  one  point  per  hop  at  the  mid-path.  It  was 
felt  that  improvements  could  be  made  by  calculating  the  absorption  on  the  up  leg  and  the  down 
leg  of  each  hop,  separately  for  each  ray.  This  involved  a  relatively  simple  calculation  of  the 
location  of  the  point  where  each  ray  penetrated  the  90  km  altitude.  Then  the  geographic 
coordinates  of  those  points  were  computed  and  the  contribution  to  the  total  absorption  was 
calculated  by  adding  the  absorption  at  each  passage  point  through  the  D-region.  Depending  on 
the  time  of  day  and  ray  path  direction,  this  could  amoimt  to  several  decibels. 

The  interface  between  PRISM  and  SMART  has  worked  very  well.  With  the 
modifications  to  the  absorption  routine  and  the  addition  of  the  ITU  noise  program,  propagation 
performance  maps  were  generated.  The  first  maps  were  centered  on  Hanscom  AFB,  using  the 
location  of  that  sounder  as  the  center  of  the  map.  The  output  of  this  program  was  compared  to 
the  full  HASEL  ray  tracing  and  the  results  compared  very  well.  Routinely,  two  maps  were 
available,  one  for  received  signal  level  and  the  other  for  signal-to-noise  ratio. 
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These  ray  tracing  calculations  were  initially  carried  out  in  one-degree  steps  in  elevation 
going  from  1°  to  90°  or  imtil  the  ray  penetrated  the  ionosphere,  and  azimuth  in  one-degree  steps, 
typically  from  0°  to  360°.  The  run  time  was  somewhat  long  and  it  was  decided  to  increase  the 
azimuthal  step  size  to  2°.  This  shortened  the  run  time  by  a  factor  of  two,  but  because  of  the 
spreading  of  fee  rays,  some  areas  at  fee  longer  ranges  were  left  unfilled. 

At  this  point  it  was  decided  to  extend  these  ray  tracings  out  to  fee  full  two-hop  coverage. 
To  do  this  fee  absorption  for  fee  two-hop  situation  ^  to  be  modified.  Wife  these  very  long 
ranges  many  of  fee  coverage  cells  (initially  1°  latitude  by  1°  longitude)  were  left  empty,  where 
no  ray  fell  on  these  cells  and  resulted  in  a  “Swiss  cheese”  appearance  on  fee  coverage  maps. 
Increasing  fee  PRISM  cell  size  to  2'  by  2°  mitigated  this  problem.  This  also  shortened  fee 
processing  time  and  now  many  more  of  fee  previously  empty  cells  had  rays  in  them.  Another 
technique  was  developed  to  eliminate  these  stray  holes  in  fee  propagation  pattern.  A  simple 
algorithm  was  used  to  determine  if  a  particular  empty  cell  had  filled  neighbors.  If  fee  answer 
was  yes,  fee  empty  cells  were  filled,  using  fee  average  of  fee  surrounding  cells.  Using  these  two 
techniques  together  produced  quite  acceptably  looking  maps.  Finally,  a  three-color  scheme  was 
developed  to  indicate  good,  average  and  poor  propagation  conditions.  This  was  intended  as  a 
guide  for  fee  unsophisticated  user. 

2.3  DIGISONDE  DATABASE 

2  J.l  Archiving  the  Digisonde  Database 

The  process  of  archiving  of  Millstone  Hill  data  continued.  Data  covering  fee  following 
time  periods  was  archived  on  4  CD-ROMs  entitled: 

98280  to  98313  (October  07, 1998  to  November  09, 1998) 

98314  to  98328  (November  10, 1998  to  November  24,  1998) 

98329  to  98360  (November  25, 1998  to  December  26, 1998) 

98361  to  99073  (December  27, 1998  to  March  14, 1999). 

From  2100  UT  on  98316  (November  12,  1998)  to  1815  UT  on  98317  (November  13, 
1998)  and  also  from  1600  UT  on  98332  (November  28,  1998)  to  0200  UT  on  98333  (November 
29,  1998),  fee  data  includes  a  mixture  of  Millstone  Hill  and  Hanscom  ionograms  with 
overlapping  measurement  times.  This  resulted  because  of  a  campaign  that  needed  data  from 
both  stations.  The  Station  ID  in  fee  preface  of  fee  Hanscom  data  was  changed  from  242  to  042 
(fee  Millstone  Hill  Station  ID)  wife  software  written  for  that  purpose,  to  avoid  repeated  error 
messages  when  fee  data  was  edited  wife  ADEP.  As  usual  before  archiving  Digisonde  data  on 
CD-ROM,  all  fee  data  were  run  through  ADEP  to  detect  corrupted  data  feat  ADEP  cannot 
handle.  All  corrupted  data  were  removed. 

The  Sondrestrom  drift  data  for  fee  time  periods  98201  to  98352  (July  20,  1998  to 
December  18,1998)  and  98353  to  99104  (December  19,  1998  to  April  14,  1999)  were  also 
archived  on  two  CD-ROMs.  A  few  minor  errors  were  foimd  in  fee  drift  data.  The  prefaces  for 
fee  data  in  file  8234J02A.dft  indicated  feat  fee  data  were  ionogram  data  instead  of  (Mft  data,  so 


10 


this  file  was  excluded  from  the  one-day  files  on  the  CD-ROMs.  This  error  and  others  were  found 
interspersed  with  the  good  data  and  were  corrected  or  removed. 

The  Millstone  Hill  Di^sonde  data  and  the  Sondrestrom  drift  data  were  continuously  sent 
to  the  UMLCAR  data  center  via  FTP  and  these  data  were  archived  on  CD-ROMs  when  sufficient 
data  had  accumulated. 

Working  back  in  time,  Millstone  Hill  data  were  archived  for  the  following  dates  in  6  CD- 

ROMs: 


98118  to  98146  (April  28, 1998  to  May  26, 1998) 

98147  to  98157  (May  27, 1998  to  June  06, 1998) 

98158  to  98176  (June  07, 1998  to  June  25, 1998) 

98 1 77  to  98208  (June  26, 1 998  to  July  27, 1 998) 

98209  to  98260  (July  28, 1998  to  September  17,  1998) 

98261  to  98279  (September  18,  1998  to  October  06, 1998). 

Also,  the  Millstone  Hill  Digisonde  data  and  the  Sondrestrom  drift  data  continue  to  be 
processed  and  archived.  Sondrestrom  drift  data  for  the  time  period  98089  to  98200  (March  30, 
1998  to  July  19,  1998)  were  also  archived  on  a  CD-ROM. 

2.3.2  Developments  of  Internet  Access  of  SAO  Data  from  Digisonde 
Stations 

The  Digisonde  online  archive  data  retrieval  package  consists  of  an  applet  providing  a 
secure  front  end  GUI  for  the  client,  i.e.,  Digisonde  Online  and  a  CGI  program  providing  on  the 
back  end,  the  search,  reading,  and  HTML  page  formatting  of  the  requested  data.  The  Digisonde 
Online  (Digisonde  OL)  was  designed  to  read  only  Standard  Archive  Output  “SAO”  data  sets,  as 
specified  during  the  Digisonde  Users  Training  Seminar  in  cooperation  with  the  World  Data 
Centers.  The  Internet-enabled  Digisonde  stations  maintain  a  local  database  of  ionospheric 
characteristics.  The  database  is  structured  into  one-day  files  in  a  Standard  Archival  Output 
(SAO)  format,  and  the  WWW  server  of  the  Digisondes  provides  software  tools  for  online 
cataloging  and  visualization  of  time  series  using  the  SAO  data  sets  in  the  database.  The 
Digisonde_OL.class  applet  was  completed. 

The  implementation  of  the  Text  and  Settings  folders  was  also  completed.  A  brief  description  of 
these  folders  follows: 
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Figure  6.  Text  folder  example 


The  text  rendering  of  the  selected  data  sets  (Figure  6)  can  be  observed  by  clicking  on  the 
tab  labeled  Text  of  die  folder  selection.  This  folder  contains  the  ASCII  text  values  of  the 
parameters  selected  and  plotted  on  the  Graphics  folder.  The  time  and  parameter  identifier  is 
given  on  the  first  line.  The  fields  are  populated  as  a  new  selection  is  made.  The  scroll  bar  allows 
the  client  to  scroll  down  to  the  end  of  the  data  set.  The  column  labeled  Processing  indicates  if  the 
selected  parameter  has  been  scaled  automatically  using  ARTIST  or  manxially  using  ADEP.  The 
text  display  is  useful  for  clients  as  a  means  to  quickly  get  at  the  actual  values  for  further  data 
manipulation.  The  text  can  be  copied  to  a  client’s  favorite  spread-sheet  and/or  data  file.  To  do 
this,  the  required  text  is  highlighted;  copied  to  an  internal  buffer,  then  the  buffer  dumped  to  the 
application  or  file. 
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New  data  set 


Figure  7.  Settings  folder  example 


The  Settings  folder  (Figure  7),  selected  by  clicking  on  the  Settings  tab,  allows  the  client 
to  change  the  color-coding  of  the  original  graphics  display.  Clicking  the  left  button  on  the  mouse 
and  running  the  mouse  cursor  over  the  horizontal  an^or  vertical  color  bars  changes  the  item’s 
color.  The  Graphics  item  is  selected  using  the  pull  down  option  above  the  color  bars.  The 
changes  are  submitted  by  pressing  the  “Set  color.”  Pressing  the  “Default  Color”  button  is  used 
for  the  original  default  color  settings.  This  folder  also  allows  the  client  to  choose  a  new  data  set. 
Simply  input  the  required  year  and  day  number  and  press  “New  data  set”  button  to  submit  the 
request.  The  new  data  set  is  loaded  automatically  and  die  Text  and  Graphics  folders  are  reset. 

While  the  applet  has  internal  settings  that  determine  security  permissions,  the  actual 
delegation  or  setup  of  the  security  model  needs  to  be  revisited  by  the  system  administrators. 

As  a  detailed  overview  on  how  this  system  works,  the  Digisonde  OL  applet  can  be 
downloaded  from  the  web  as  part  of  an  HTML  document.  It  begins  its  life  on  the  web  page  as  a 
single  button  labeled  “DigisondeDB”.  Once  this  button  is  pressed  the  applet  calls  a  CGI  program 
(DBscan)  to  find  the  most  recent  data  set  available  to  the  client.  It  then  unpacks  the  data  and 
formats  a  return  response  within  an  HTML  document.  At  this  point  the  client  would  be  presented 
with  the  main  GUI  frame.  This  main  frame  consists  of  several  components.  These  are: 
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A  station  name,  data  set  year,  and  day  number  header  at  the  top  of  the 
1^1  ^plet  screen.  This  allows  for  rendering  multiple  views  of  several  data 
"■  sets  from  the  same  station  and/or  views  of  data  sets  from  other 

stations.  Every  time  the  ThgisondeDB”  button  on  the  originai  HTML 
page  is  pressed  a  new  independent  frame  would  be  created. 


A  scrollable  parameter  selection  text  field  will  be  supplied  to  display 
die  parameters  available  for  plotting  from  the  current  data  set  Double 
clicking  on  any  parameter  loads  that  ^dfic  data  set  into  die  text 
folder  and  plots  it  on  the  graph  folder.  A  limit  of  five  parameters  will 
be  allowed  on  the  same  graph.  An  error  message  would  be  given  if  an 
attenqit  is  made  to  load  any  more  than  five  parameters. 


A  set  of  control  buttons  will  be  supplied.  These  buttons  will  consist  of. 

Clear:  This  button  allows  the  client  to  remove  selected  parameters 
individually  or  all  die  selected  parameters. 

Download:  This  button  initiates  a  download  using  options  in  the 
originai  HTML  page  from  which  the  applet  originated. 

Help:  As  the  name  suggests  this  opais  a  help  box  with  information  on 
how  to  use  the  Digisonde  OL  front  end. 

Quit:  Closes  the  fiame  and  releases  systems  resources. 
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Figure  8.  Graph  folder  example 

The  Graph  folder  (Figure  8)  is  the  default  folder  of  three  folders  that  can  be  selected  from 
the  Graphics,  Text,  and  Settings  tabs.  When  a  new  frame  is  generated  the  graphics  view  of  the 
data  set  is  always  loaded  by  default.  This  view  plots  the  selected  parameter  values  over  a  24-hour 
period.  Universal  time  is  represented  along  the  X-axis,  while  on  the  Y-axis  ,  the  frequency  scale 
(MHz)  is  on  the  left  and  height  scale  (km)  is  on  the  right. 

Below  the  plot  is  the  text  representation  of  the  selected  parameters  for  specific  times. 
When  the  mouse  cursor  is  moved  over  the  plot,  the  corresponding  time  and  parameter  values  for 
the  cursor  position  will  be  updated  in  the  text  fields.  When  no  parameters  are  selected  the  main 
identifiers  in  the  text  area  are  labeled  Time,  Graph  1,  Graph  2,  Graph  3,  Graph  4  and  Graph  5. 
As  each  parameter  is  selected  the  labels  are  updated.  As  an  example,  when  foF2  is  selected  as  the 
first  parameter,  the  label  Graph  1  will  be  changed  to  foF2,  and  the  adjacent  text  field  will  be 
active.  Parameter  values  in  these  text  fields  share  the  same  coloring  code  as  the  plotted  curves 
for  easy  identification  between  graph  and  text. 

These  items  have  been  successfully  implemented  in  Java.  The  Digisonde  OL  applet  used 
the  Java.net  and  Sim.net  class  libraries  to  provide  the  Application  Programming  Interface  (API) 
necessary  to  access  Internet  and  Web  resources.  Hence  these  APIs  provide  a  similar  approach  to 
distributed  computing  as  a  Web  Browser  in  that  the  applet  is  able  to  use  the  Common  Gateway 
Interface  (CGI)  to  invoke  the  program  “Dbscan.exe”  on  an  HTTP  server.  The  Dbscan.exe 
currently  accepts  two  arguments,  the  year  and  day  number.  It  does  a  search  for  the  requested  data 
set,  unpacks  the  appropriate  SAG  file,  formats/packs  the  data  into  an  HTML  page  and  sends  the 
path/filename/and  data,  back  to  the  client.  The  DBscan.exe  program  is  written  in  C-H-  and  is 
maintained  on  the  Web  server  generally  in  the  cgi-bin  directory. 
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A  local  database  of  ionospheric  characteristics  was  designed  to  allow  for  direct  access  to 
the  time  series  over  the  Internet  An  archiving  scheme  was  developed  to  provide  and  manage  the 
required  daily  collections  of  ionogram  characteristics.  Catalog  browsing  software  was  developed 
to  dynamically  create  availability  catalogs  for  the  database  contents  and  allow  for  data  download. 
A  scheme  was  developed  to  ensure  restricted  access  to  the  secure  data  in  the  database.  Client- 
side  software  was  developed  to  provide  for  interactive  visualization  of  the  ionospheric 
characteristics  and  data  flow  control.  Data  exchange  protocols  between  client  and  server  were 
modified  to  provide  for  a  greater  number  of  characteristics  for  display,  proper  station 
identification,  and  an  indicator  of  data  quality  validation.  A  scheme  was  developed  to 
automatically  search  for  the  latest  available  data  in  the  archive  on  the  first  access. 

Finally,  the  introduction  of  the  newer  SAO  4.2  format  required  ftmctional  testing  of  the 
local  Digisonde  SAO  database  support  software.  The  display  software  had  to  be  modified  to 
allow  proper  selection  and  visualization  of  the  new  characteristics  that  were  introduced  in  SAO 
4.2. 

The  characteristic  names  on  the  display  were  brought  into  agreement  with  URSI 
recommendations.  Also,  these  tests  showed  an  inconsistency  of  operational  controls  on  the 
display  window,  which  was  remedied  by  introduction  of  a  callback  routine  monitoring  the 
window  “close”  event. 

2.4  DIGISONDE  NETWORK  SUPPORT 
2.4.1  ARTIST4  Upgrades 

For  the  planned  hardware  shipment  to  Chile,  techniques  for  archiving  the  16-channel 
“raw”  data  from  the  Digisonde  256  Processor  Output  computer  were  investigated.  The  solution 
arrived  at  was  to  modify  the  FEND.EXE  program  in  the  ARTISTS  software  to  produce  an 
additional  file  type,  .16c  files,  when  PI  =  6  in  the  ionogram  program.  The  .16c  files  would  be 
sent  to  the  ARTIST4  computer  where  Dispatcher  would  take  care  of  archiving  them. 

Specifications  for  ARTIST4  Phase  2  hardware  were  developed.  The  major  changes  were 
the  substitution  of  an  IOMEGA  JAZ  drive  in  place  of  the  Tandberg  TDC  3820  tape  drive,  and 
the  addition  of  a  B&B  watchdog  timer  card.  Also  an  additional  32  megab34es  of  memory  were 
added  to  make  a  total  of  64  Megabytes.  Additional  ARTIST4  computers  meeting  the  above 
specifications  were  ordered  along  wiA  memory,  watchdog  timer  cards,  and  JAZ  drives  to  update 
the  earlier  five  computers. 

To  augment  this  effort  a  diagnostic  procedures  document  for  the  ARTIST4/ARTIST3/ 
network  coimection  was  developed.  The  ARTISTS  portion  of  Web-based  remote  control 
software  was  developed. 

Y2K  fixes  were  made  to  the  processor  chassis  firmware,  some  of  the  ARTISTS 
programs,  and  the  ARTIST4  program.  Y2K  improvements  were  made  to  the  Timex  module, 
which  is  used  in  the  DARTIST3.Exe  and  Dps  Control  programs.  A  major  activity  was  to  ensure 
that  all  of  the  USAF  DISS  installations  would  be  Y2K  compliant  before  the  end  of  the  1998. 
UMLCAR  personnel  participated  in  the  official  USAF  Y2K  test  for  the  ARTIST4  upgraded 
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DISS.  The  test  was  performed  on  28  October  1998  at  the  Wallops  Island,  VA  DISS  site  imder 
the  supervision  of  MSgt.  Robert  Tadena,  AFMC  SMC  Det  1  l/CBD.  The  tests  were  successful. 

Much  Web  software  was  developed  for  purposes  of  control  of  the  DISS  over  the  Web 
and  more  user-friendly  display  of  data.  A  Netscape  FastTrack  server  was  added  to  the  ARTIST4 
package  in  place  of  the  Microsoft  Windows  Peer  Web  Service  Web  Server.  Changes  were  made 
to  the  ARTIST4  cgi  scripts  to  be  compatible  with  the  Netscape  program. 

Also  a  non-linear  frequency  scan  was  developed.  Primarily,  this  required  changes  to  the 
processor  chassis  card  2  firmware  to  select  frequencies  based  on  a  table.  Much  thought  went 
into  the  non-linear  frequency  table  by  several  individuals  based  on  bandwidth  of  the  Digisonde 
and  critical  frequency  resolution  requirements.  The  ARTISTS  program  that  inputs  data  from  the 
processor  chassis,  FEND.EXE,  was  modified  to  accept  the  data  and  pass  it  on  to  DARTIST3  and 
ARTIST4.  The  Web  portion  of  the  Web-based  remote  control  software  was  developed. 

A  full  security  scheme  was  devised  and  implemented.  The  scheme  involved 
configuration  of  Netscape  FastTrack  Server,  Windows  NT,  Microsoft  Peer  Web  Services,  and 
Microsoft  Remote  Access  Service  (RAS).  Different  levels  of  security  were  provided  for  Web 
access,  FTP  access,  and  local  access  depending  on  the  user’s  permissions. 

UMLCAR  personnel  went  to  SMALC  in  Sacramento,  CA  on  16  November  98  for  the 
purpose  of  training  Dan  Rivas  and  Clay  Dukes,  SMALC  in  Sacramento,  on  the  Phase  n  of  the 
ARTIST4  upgrade  to  the  DISS.  Further  training  and  validation  of  procedures  and  documentation 
occurred  on  30  November  98  at  the  Eglin  AFB,  FL,  DISS  site.  Dan  Rivas  participated  in  the 
upgrade  of  Phase  n  of  the  ARTIST4.  Similar  training  and  validation  occurred  on  2  December  98 
at  the  Ramey  Solar  Observatory,  PR,  DISS  site.  Dan  Rivas  and  Bruce  Andrews,  from 
Tobyharma,  PA,  Army  Depot  participated  in  a  complete  ARTIST4  installation.  On  27  December 
98,  UMLCAR  persormel  went  to  the  Bermuda  DISS  site  and  did  a  complete  ARTIST4 
installation.  Discussions  were  held  with  the  remaining  NASA  persormel  concerning  the  changes, 
which  would  have  to  be  made  to  accommodate  the  reduced  NASA  presence  on  the  island. 
Instructions  were  prepared  and  support  was  provided  to  the  local  persormel  at  Goose  Bay, 
Canada,  and  Narssarssuaq,  Greenland,  so  that  the  local  persormel  could  perform  ARTIST4  Phase 
n  upgrades  to  their  DISS  installations  that  had  previously  been  upgraded  by  UMLCAR 
persormel  to  ARTIST4  in  August  98.  Similarly,  instructions  and  support  were  provided  for  local 
persormel  at  Sondrestrom  and  Qaanaaq,  Greenland,  so  that  the  local  persormel  could  perform 
complete  ARTIST4  installations.  Complete  ARTIST4  kits  were  also  prepared  and  shipped  to 
Korea,  Australia,  England,  Italy,  Alaska,  Sondrestrom,  Qaanaaq,  Puerto  Rico,  and  Bermuda. 

2.4.2  DISS  Network  Support 

A  system  for  remote  control,  status  visualization  and  diagnostics  of  DISS  over  the 
standard  HTTP  layer  was  designed.  A  study  of  dynamic  control  mechanisms  regarding  their 
conummications  overhead,  browser  compatibility,  hard  disk  activity  and  response  time  was 
conducted.  The  control  schemes  "server  push"  and  "client  pull"  were  evaluated.  Both  schemes 
were  used  to  develop  alternative  software  solutions  and  they  were  tested  side-by-side.  The  server 
push  technology  proved  to  be  more  efficient  in  reporting  slow,  uneven,  or  asynchronous  data 
updating.  However,  imbridgeable  compatibility  problems  were  met  because  the  browser 
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software  companies  attempt  to  undermine  each  other's  development  by  enforcing  incompatible 
interfaces.  For  this  reason  the  client-pull  solution  was  selected  for  the  final  development.  To 
exclude  the  server  fi'om  extensive  transaction  logging  activity  in  the  client-pull  scenario,  the 
Microsoft  server  software  was  changed  to  Netscape’s  FastTrack  server  with  a  wider  range  of 
features.  The  developed  software  provides  self-updated  displays  of  system  status  and  near  real¬ 
time  ARTISTS  screen  capture,  as  well  as  the  means  for  remote  commanding  of  the  DISS. 

A  local  long-term  archive  of  daily  ionospheric  data  was  designed.  A  data  quality  control 
scheme  was  implemented  to  ensure  valid  data  formats.  Housekeeping  procedures  were  adjusted 
to  include  new  types  of  data  as  well  as  to  maintain  an  archive  of  data  with  format  errors.  A  new 
archival  media  was  introduced  for  easier  access  to  the  long-term  archive.  System  control 
software  was  modified  to  collect  and  display  the  media  status. 

The  new  preface  data  format  was  designed  to  accommodate  the  year  2000  transition. 
Data  handling  routines  were  modified  to  include  new  types  of  automatically  scaled  information 
and  to  operate  with  a  new  Y2K-compliant  format  of  system  settings.  Additions  to  the  DISS 
control  panel  were  made  to  include  online  ARMENU  editor  and  extended  screen  capture  history 
displays. 

Provisions  were  made  in  the  Sondrestrom  ARTIST4  configuration  so  that  local  personnel 
could  have  access  to  and  be  able  to  control  the  DISS.  Also,  provisions  were  made  so  that  several 
scientific  organizations  could  access  real-time  data. 

2.5  Digital  Sounder  for  Topside  Observations 

Special  Topic  -  Staggered  Pulse  Sequences,  Basic  Relationships 

Introduction 

The  Staggered  Pulse  Sequence  (SPS)  waveform  is  investigated  for  possible  deployment 
in  topside  soimders.  It  consists  of  pseudo-randomly  spaced  single  pulses,  each  with  a  pulse 
width: 


W  =  1  /  (receiver  bandwidth). 

The  pulse  width  should  be  sufficiently  narrow  to  resolve  multiple  echo  sources  separated  by  a 
resolution  distance: 


AR  =  W*c/2, 

where  c  is  the  speed  of  light.  The  pseudo-random  spacing  is  such  that  echoes  from  previously 
transmitted  pulses  are  received  between  the  transmission  of  subsequent  pulses.  Since  echoes 
may  come  from  any  range,  some  of  the  pulse  echoes  are  obscined  (in  a  monostatic  configuration) 
by  the  subsequent  transmitter  pulses,  so  the  pulse  pattern  has  to  be  such  that  an  equal  number  of 
echoes  are  available  from  all  ranges.  Another  complication  of  this  scheme  is  that  at  times 
corresponding  to  an  echo  from  one  source  from  pulse  “i”,  another  echo  from  a  second  source, 
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from  a  different  range  and  in  response  to  pulse,  “j”,  may  be  arriving  simultaneously.  Therefore 
the  optimum  pattern  for  these  pulses  would  provide  a  nearly  equal  number  of  clear  echoes  from 
every  range  and  will  create  a  minimum  number  (or  assuming  all  ranges  are  of  equal  importance, 
an  evenly  distributed  number)  of  spurious  echoes  when  sources  exist  at  more  than  one  range. 

The  signal  processing  gain  is  realized  by  coherently  integrating  the  received  echoes. 
Since  the  moving  sources  results  in  a  constantly  changing  phase  in  the  echo  signals,  this 
integration  must  be  performed  in  the  frequency  domain.  Since  the  pulse  echoes  arrive  at  variable 
(pseudo-random)  intervals,  the  spectrum-integration  task  becomes  one  of  estimating  the 
spectrum  at  each  range  from  a  set  of  non-imiformly  spaced  samples. 

There  are  many  parameters,  each  of  which  can  be  specified  over  a  wide  range,  and  the 
problem  of  optimizing  the  design  must  be  thought  out  before  attempting  to  define  particular 
parameters.  The  design  must  include: 

1 .  The  determination  of  number  of  pulses  in  the  sequence. 

2.  The  duty  factor  for  those  pulses. 

3.  The  desired  minimum  number  of  clear  echoes  from  each  range  to  be  integrated,  thereby 
realizing  a  particular  signal  processing  gain. 

4.  The  length  of  the  sequence,  which  gives  the  time  over  which  echo  pulses  can  be 
integrated. 

Spectral  Characteristics 


Of  course,  another  significant  part  of  the  design  was  to  specify  the  spectrum-integration 
technique,  which  will  be  addressed  later.  Suffice  it  to  say  that  one  of  the  requirements  coming 
from  the  spectrum-integration  task  is  that  there  be  a  minimum  gap  between  available  samples  at 
every  range  and  that  the  spacing  of  the  available  samples  is  as  random  as  possible,  i.e.,  all  auto 
correlation  lags  are  equally  well  represented.  Further  considerations  are  that  the  pulses  need  to 
span  a  sufficient  time  interval  (the  coherent  integration  time  or  CIT)  to  provide  the  desired 
Doppler  resolution  when  the  N  pulses  are  coherently  integrated. 

Basic  Relationships 

The  number  of  clear  echoes  (C)  is  proportional  to  the  number  of  echoes  transmitted  (N) 
and  the  duty  factor  of  the  transmissions  (D).  If  the  total  time  of  the  pulse  sequence  is  M,  in  units 
of  time  slots,  each  of  width  W,  from  the  range  resolution  considerations  described  above,  the 
duty  factor  can  be  expressed  as: 


D  =  N/M. 
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Since  the  reception  of  X  =  DN  echo  pulses  will  be  blocked  due  to  subsequent 
transmissions,  the  expected  number  of  clear  pulse  echoes  at  each  range  is  approximately: 

C  =  (N-DN)  =  N(1-D). 

Next,  the  number  of  spurious  echoes  from  other  ranges  that  contaminate  samples  made  at 
times  corresponding  to  the  delayed  pulse  pattern  received  from  range  R  was  considered.  In 
general,  if  there  are  many  echo  sources  at  different  ranges,  echoes  from  other  ranges  only  occur 
infrequently  (based  on  the  duty  factor,  D)  at  sample  times  corresponding  to  the  range  being 
processed.  First,  consider  the  case  of  two  distinct  echo  sources,  &e  wanted  echo  at  range  R 
(assuming  we  are  currently  echo  detecting  for  range  R)  and  a  spurious  source  at  range  S.  If  there 
was  more  than  one  echo  from  range  S  they  integrate  coherently  (although,  in  general  at  a 
different  spectral  line)  when  a  full  spectral  integration  is  made.  The  probability  of  an  echo  from 
range,  S,  falling  at  a  time  interval  corresponding  to  another  echo  from  range  R  is  equal  to  the 
duty  factor,  D.  That  is,  the  lower  the  density  of  pulses  traveling  to  and  from  the  echo  sources, 
the  lower  the  joint  probability  of  two  echoes  arriving  simultaneously.  However,  since  we  are 
discussing  the  particular  time  at  which  an  echo  from  R  is  being  received,  it  is  not  the  joint 
probability,  but  the  conditional  probability  that  detennines  how  likely  it  is  that  an  echo  from  S  is 
also  present.  As  was  shown  above,  there  were  N(l-D)  clear  echoes  from  R,  and  in  them  will  be 
N(1-D)D  spurious  echoes  from  range  S.  Therefore,  the  ratio  of  wanted  echoes  to  spurious 
echoes  from  any  other  range  is  simply  D.  Therefore,  it  can  be  conclude  that: 

1 .  If  a  processing  gain  G  is  desired,  C  pulses  must  be  received  to  achieve  this. 

2.  To  reduce  spurious  interference  the  time  M  over  which  the  C  pulses  are  received 
should  be  stretched  out  as  much  as  possible  to  lower  the  duty  factor. 

Assuming  that  the  Doppler  analysis  range  is  very  large  compared  to  actual  signal 
Doppler,  it  would  be  advantageous  to  apply  an  interpulse  phase  code  at  transmission,  so  that  the 
leakage  observed  at  one  range  due  to  pulses  from  another  changes  phase  rapidly.  The  phase  of 
the  wanted  echoes  must  be  inverted  upon  reception  to  compensate  for  Ae  phase  reversals 
modulated  onto  the  transmitted  signal;  then  the  cross-talk  echoes  will  arrive  with  random  phase. 
A  near  optimmn  pseudo-random  code  has  been  generated  for  the  768  time-slot  code  that 
improves  the  cross-talk  leakage  by  randomizing  the  spectrum.  This  improves  the  narrowband 
spectral  signal  (from  the  wanted  range)  to  leakage  ratio. 

The  practical  limits  on  low  duty  factor  are  the  desired  speed  of  the  measurement  and  the 
temporal  coherence  of  the  medium  propagating  the  signal.  As  noted  above,  in  the  section  on 
“Spectral  Characteristics,”  the  Doppler  resolution  and  Doppler  range  (respectively,  the  inverse  of 
the  err  and  the  inverse  of  the  closest  pulse  spacing)  are  affected  by  the  length  of  the  sequence 
(M).  If  a  very  high  Doppler  range  is  required  it  may  be  necessary  to  allow  pulses  to  occupy 
adjacent  time  slots  in  groups  of  1,  2,  3  or  possibly  a  greater  number  of  transmitter  pulses. 
However,  here  only  the  specific  case  of  single  pulse,  spaced  by  at  least  two  receiver  (quiet)  time 
slots  was  considered.  Two  time  slots  are  needed  since  there  is  generally  a  receiver  recovery  time 
of  at  least  one  time  slot  before  echoes  can  be  received. 
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Next,  the  possibility  of  many  spurious  echoes  arriving  from  different  ranges,  Si,  Sa-.-.Sp 
was  considered.  In  general  the  phase  and  its  first  derivative,  the  Doppler  shift,  of  spurious 
echoes  amving  from  various  ranges  are  independent  and  therefore  sum  up  incoherently  in  the 
spectrum-integration  process,  so  that  the  effect  of  multiple  sources  is  only  large  if  too  many 
ranges  provide  echoes.  Now  it  is  necessary  to  compare  the  range  leakage  characteristics  of  SPS 
to  conventional  pulse/Doppler  radars  transmitting  a  pulse  train  (e.g.,  64  or  128  pulses)  and 
coherently  integrating  the  received  echoes.  In  most  radar  geometries  there  is  a  maximum  range 
(Rmax)  that  provides  echoes,  or  beyond  which  the  echoes  are  too  weak  to  be  detected.  Usually 
the  pulse  repetition  rate  is  chosen  such  that  all  possible  pulse  echoes  are  received  before  another 
pulse  in  the  pulse  train  is  transmitted,  thus  avoiding  range  ambiguity  and  range  aliasing.  A  ratio 
of  M  to  Rmax  was  quantified  that  describes  how  the  pulse  density  (the  duty  factor)  can  be  varied 
all  the  way  to  the  case  of  the  conventional  radar  as  a  limiting  case  of  the  SPS  parameters.  The 
conventional  radar  sends  one  pulse  per  Tmax,  the  time  delay  associated  with  the  echo  from 
Rmax>  where  Tmax  -  Rmax  /(c/2). 


In  the  case  of  SPS,  the  transmitted  code  sequence  is  much  longer  than  the  number  of 
ranges  of  interest  (M(c/2)  >  Rmax)-  Therefore  as  echoes  are  received  near  the  end  of  the  SPS 
sequence,  it  is  not  likely  to  receive  any  very  long  delayed  echoes  from  the  pulses  transmitted 
near  the  beginning  of  the  sequence.  In  the  analysis  shown  below,  the  search  for  spurious  echoes 
to  ranges  less  than  256  have  been  cut  off.  This  tradeoff  has  already  been  taken  into  account  in 
the  expressions  for  spurious  echo  density  above  in  that  the  lower  duty  factor  results  in  a  longer 
sequence  and  the  longer  sequence  has  a  lower  density  of  sptuious  echoes  at  each  range. 
Therefore  of  the  number  of  previously  transmitted  pulses  (iD,  where  “i”  is  index  of  the  time  slot 
at  which  a  clear  echo  is  received),  only  a  certain  number  of  them  (on  the  average  over  the  entire 
sequence,  Rmax  /M)  *  N  )  may  end  up  as  spurious  echoes.  However,  if  the  ratio  of  M  to  Rmax  is 
high  a  large  number  of  spurious  echoes  will  never  be  seen.  Of  course  one  of  the  limiting  cases 
of  a  very  low  duty  factor  (which  results  in  a  very  large  ratio  of  M/Rmax)  is  the  conventional  low 
PRF  design  where  spurious  echoes  are  never  seen  because  they  come  from  ranges  that  do  not 
exist  or  caimot  be  detected. 

3.  SUMMARY 

Progress  has  been  made  on  modeling  the  equatorial  airglow  using  a  relatively  simple  set 
of  differential  equations.  This  approach  has  made  it  possible  to  generate  electron  density  profiles 
and  emission  profiles  as  a  function  of  time,  particularly  after  simset.  These  calculations  generate 
the  background  electron  density  and  630.0  nm  emission  as  a  function  of  time  in  which  the 
depletions  are  embedded.  The  equations  controlling  the  electron  density  variation  with  height 
and  time  require  a  vertical  velocity  profile  input.  This  means  that  by  comparing  the  calculated 
electron  density  profiles  with  the  measured  ones,  it  is  possible  to  generate  a  reasonable  velocity 
profile.  These  results  ultimately  can  be  used  to  determine  the  associated  electric  field  in  the  F- 
region  that  drives  the  layer  up  or  down. 

A  new  project  was  added  to  this  year’s  effort  involving  propagation  analysis  over  a  wide 
coverage  area.  Project  SEND  was  designed  to  evaluate,  using  the  best  updated  ionospheric 
models  and  rapid  ray  tracing  software,  the  effectiveness  of  a  particular  communications  system 
in  providing  service  within  the  coverage  area.  The  program  allows  users  to  effectively  manage 
frequency  and  power  of  their  systems  in  an  almost  real-time  environment.  Considerable  effort 
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was  required  to  evaluate  different  ray  tracing  schemes  and  ionospheric  modeling  techniques. 
Additions  to  the  software  models  were  made  to  achieve  the  desired  goals.  It  only  remains  to  test 
the  operational  software  and  evaluate  its  performance  in  the  real  world.  On  the  basis  of  these 
tests,  some  refinements  to  the  code  are  expected. 

Considerable  effort  was  expended  this  year  in  getting  the  Air  Force’s  DISS  network  “on 
line,”  addressing  questions  on  data  format,  security  and  ease  of  access.  The  Y2k  issue  was 
addressed  and  software  modifications  made  and  tested  to  demonstrate  that  the  Digisonde 
network  would  operate  successfully  through  the  year  2000  transition.  The  ARTIST4  program 
was  completed  and  installed  in  several  of  the  DISS  systems  in  the  US  and  elsewhere.  TTie  SAO 
formatting  and  the  ARTIST4  package  make  the  Digisonde  network  very  effective  in  providing  a 
global  sounding  system  that  serves  both  the  US  Air  Force  and  the  scientific  researchers  who 
require  these  data. 

Finally,  as  part  of  die  design  of  the  new  topside  sounding  system  a  new  waveform  was 
introduced  involving  staggered  pulse  sequencing.  The  effectiveness  of  this  proposed  waveform 
was  demonstrated  in  terms  of  its  range  and  Doppler  performance. 
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